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We show that in theories of gravity that add quadratic curvature invariants to the Einstein- 
Hilbert action there exist expanding vacuum cosmologies with positive cosmological constant which 
do not approach the de Sitter universe. Exact solutions are found which inflate anisotropically. This 
behaviour is driven by the Ricci curvature invariant and has no counterpart in the general relativistic 
limit. These examples show that the cosmic no-hair theorem does not hold in these higher-order 
extensions of general relativity and raises new questions about the ubiquity of inflation in the very 
early universe and the thermodynamics of gravitational fields. 
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I. INTRODUCTION 

The inflationary universe is the central cosmological 
paradigm which astronomical observations aim to test, 
and by which we seek to understand how the universe 
might have evolved from a general initial condition into 
its present state of large-scale isotropy and homogeneity 
together with an almost flat spectrum of near-Gaussian 
fluctuations. The essential feature of this inflationary 
picture is a period of accelerated expansion during the 
early stages of the universe Q] . The simplest physically- 
motivated inflationary scenario drives the acceleration 
by a scalar field with a constant potential, and the lat- 
ter can also be described by adding a positive cosmo- 
logical constant to the Einstein equations. In order to 
understand the generality of this scenario it is impor- 
tant to determine whether universal acceleration and 
asymptotic approach to the de Sitter metric always oc- 
curs. A series of cosmic no-hair theorems of varying 
strengths and degrees of applicability have been proved 
to demonstrate some necessary and sufficient conditions 
for its occurrence 0, 0> LH M E| • Similar deductions 
are possible for power-law [3, |8| an d intermediate in- 
flationary behaviour, where accelerated expansion is 
driven by scalar-field potentials that have slow exponen- 
tial or power-law fall-offs, but we will confine our dis- 
cussion to the situation that occurs when there is a pos- 
itive cosmological constant, A > 0. So far, investiga- 
tions have not revealed any strong reason to doubt that, 
when A > and other matter is gravitationally attrac- 
tive, any stable, ever-expanding general-relativistic cos- 
mological model will approach isotropic de Sitter infla- 
tion exponentially rapidly within the event horizon of any 
geodesically-moving observer. Similar conclusions result 
when we consider inflation in those generalisations of gen- 
eral relativity in which the Lagrangian is a function only 
of the scalar curvature, R, of spacetime. This similarity 
is a consequence of the conformal equivalence between 
these higher-order theories in vacuum a nd g eneral rela- 
tivity in the presence of a scalar field [HJ, EH E2 • I n this 
paper we will show that when quadratic terms formed 



from the Ricci curvature scalar, R^ V R^ U are added to the 
Lagrangian of general relativity then new types of cosmo- 
logical solution arise when A > which have no counter- 
parts in general relativity. They inflate anisotropically 
and do not approach the de Sitter spacetime at large 
times. We give two new exact solutions for spatially ho- 
mogeneous anisotropic universes with A > which pos- 
sess this new behaviour. They provide counter-examples 
to the expectation that a cosmic no-hair theorem will 
continue to hold in simple higher-order extensions of gen- 
eral relativity. Other consequences of such higher-order 
theories have been studied in [l3L fLU Il5l | . The presence 
of such quadratic terms as classical or quantum correc- 
tions to the description of the gravitational field of the 
very early universe will therefore produce very different 
outcomes following expansion from general initial condi- 
tions to those usually assumed to arise from inflation. 
This adds new considerations to the ap plic ation of the 
chaotic and eternal inflationary theories |16| in conjunc- 
tion with anthropic selection |17| . 

We will consider a theory of gravity derived from an 
action quadratic in the scalar curvature and the Ricci 
tensor. More specifically, ignoring the boundary term, 
we will consider the D-dimensional gravitational action 



S G = -!- / d D x^/\g\ (R + aR 2 + 0R„ u R" v - 2A) . 
* K Jm 

(1) 

Variation of this action leads to the following generalised 
Einstein equations (see, e.g., 01): 



G^ + Q^+Ag^ = «T /1I/ , (2) 



where is the energy-momentum tensor of the ordi- 
nary matter sources, which we in this paper will assume 
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to be zero, and 



(3) 



3v = 



2aR ( R pv - jRg^ 



(2a + /3) {g^U-V^ v )R 



+I3U ( R„ v - -Rg^ 



(4) 



with □ = V M V M . The tensor <I> M „ incorporates the de- 
viation from regular Einstein gravity, and we see that 
a = (3 = implies = 0. 

First, consider an Einstein metric, so that R^ — \g pu - 
This is a solution of eq. J2J) with T M „ = provided that 

A = ±[(D-4)(Da + f3)\ + (D-2)}. (5) 

Hence, when D = 4 any Einstein space is a solution to 
eq.® provided that A = (D — 2)A/2. In particular, if 
A > 0, de Sitter spacetime is a solution to eq.©. If 
A = 0, we need A = and de Sitter spacetime cannot be 
a solution. 

Now consider solutions to eq.(J2J which are non- 
perturbative and a and (3 are not small. We know that 
solutions with (3 = 0, a ^ are conformally related to 
Einstein gravity with a scalar field <fi = ln(l + 2aR) 
that possesses a self-interactio n p otential of the form 
V{cj)) = (e* - l) 2 /4a, EH EH El, and their inflation- 
ary behaviours for small and large \<j)\, along with that 
of theories derived from actions that are arbitrary func- 
tions of R, are well understood. However, there is no 
such conformal equivalence with general relativity when 
(3 7^ and cosmologies with A > can then exhibit quite 
different behaviour. 



II. THE FLAT DE SITTER SOLUTION 

First consider the spatially-flat de Sitter universe with 
metric 



ds 2 dS = -dt 2 + e 2Ht ^~ 2 ' ^ ' ^ 



(da; 2 + dy 2 + dz 2 ) , H 



(6) 

The stability of this solution in terms of perturbations of 
the scale- factor depends on the sign of (3a + (3). In 4D, 
we can use the Weyl invariant and the Euler density, E, 
defined by |io| . 

n /iiwpa _ p -D^vpa nn r>pv lip 2 

E = R^ P aR^ pa -iR^R^ + R 2 , (7) 

to eliminate the quadratic Ricci invariant in the action, 
since 

aR 2 + f3R^R^ = ±(3a + (3)R 2 + § {C^ pa C^ pa - E) . 



Since integration over the Euler density is a topological 
invariant, the variation of E will not contribute to the 
equations of motion. The Friedmann-Robertson- Walker 
(FRW) universes are conformally flat so, for a small vari- 
ation, the invariant C AI „ p(T C A " ypcr will not contribute ei- 
ther. Hence, sufficiently close to a FRW metric only the 
R 2 term will contribute. The stability of the FRW uni- 
verse is therefore determined by the sign of (3a + (3) EH • 
One can check this explicitly using eq. © . We start with 
the metric ansatz: 



ds 2 



-dt 2 + e 2b ^ (dx 2 + dy 2 + dz 2 ) 



H 



A 

3' 




and note that in 4D the trace of eq. @ reduces to 

-i? + 2(3a + /?)QR + 4A = 0, (8) 

which can be used to determine the stability of the Ricci 
scalar. We can perturb the Ricci scalar by assuming a 
small deviation from the flat de Sitter metric of the form: 

b(t) =Ht + b ie Xlt + b 2 e X2t + e>(e 2A **), 

where b\ and bi are arbitrary constants. Eq.© implies 



(9) 



if (3a + (3) ^ 0. For (3a + (3) =0, we must have 
&i = 62 = 0. From this expression we see that if 
(3a + 0) > then the solution will asymptotically ap- 
proach the flat de Sitter spacetime as t — ► 00; however, 
for (3a + f3) < the solution is unstable. For the spe- 
cial case of J3 = 0, this result agrees with the stability 
analysis of 19]. A construction of an asymptotic series 
approximation around the de Sitter metric for the case 
13 = has also been performed [S3, EH 121 [H H3 - In the 
case of general relativity (a = (3 = 0) a number of results 
for the inhomogeneous case of small perturbations from 
isotropy and homo geneity when A > have also been 
obtained @, 0, H HM HI |H H3 . 

We see that, as long as (3a + [3) > 0, any FRW model 
sufficiently close to the flat de Sitter model will asymp- 
totically approach de Sitter spacetime and consequently 
obeys the cosmological no-hair theorem. We should em- 
phasize that only FRW perturbation modes have been 
considered here. The question of whether the flat de Sit- 
ter universe is stable against general anisotropic or large 
inhomogeneous perturbations when a / and (3 ^ is 
still unsettled. In the case of universes that are not 'close' 
to isotropic and homogeneous FRW models we shall now 
show that the cosmic no-hair theorem for A > vacuum 
cosmologies is not true: there exist ever-expanding vac- 
uum universes with A > that do not approach the de 
Sitter spacetime. 



III. EXACT ANISOTROPIC SOLUTIONS 

We now present two new classes of exact vac- 
uum anisotropic and spatially homogeneous universes of 
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Bianchi types II and Vlh with A > 0. These are new 
exact solutions of the cqns. J3J) with (a,0) ^ (0,0). 
Bianchi type II solutions: 



ds 2 = _ df 2 + e 2M 



where 



da; + -(zdy - ydz) +e bt (dy 2 +dz 2 ), 

(10) 



a 2 = 



ll + 8A(lla + 3/3) 
30(3 : 



2 _ 8A(o + 3/3) + 1 

30/3 ' 1 ' 



These solutions are spacetime homogeneous with a 5- 
dimensional isotropy group. They have a one-parameter 
family of 4-dimcnsional Lie groups, as well as an iso- 
lated one (with Lie algebras A\ 11 and A\ g , respectively, 
in Patera et al's scheme |2|j) acting transitively on the 
spacetime. An interesting feature of this family of solu- 
tions is that there is a lower bound on the cosmological 
constant, given by A min = — l/[8(a + 3/3)] = — a 2 /8 for 
which the spacetime is static. For A > A m ; n the space- 
time is inflating and shearing. The inflation does not 
result in approach to isotropy or to asymptotic evolution 
close to the de Sitter metric. Interestingly, even in the 
case of a vanishing A the universe inflates exponentially 
but anisotropically. We also note from the solutions that 
the essential term in the action causing this solution to 
exist is the (3 R^R^ '-term and the distinctive behaviour 
occurs when a = 0. The solutions have no well defined 
f3 — > limit, and do not have a general relativistic coun- 
terpart. They are non-perturbative. Similar solutions 
exist also in higher dimensions. Their existence seem to 
be related to so-called Ricci nilsolitons [30l l3l| . 

Bianchi type Vlh solutions: 



ds 2 



VI, 



= -dt 2 + dx 2 



+ e 2(rt+ax) L-2(st+ahx) dy 2 + e +2(st+ahx) ^al 



where 

„.2 



8(3s 2 + (3 + h 2 ){\ + 8Aa) + 8A/?(1 + h 2 ) 
8f3h 2 ' 



8[3s 2 + 8A(3a + (3) + 3 

zjm 2 ' 



(13) 



and r, s, a, and h are all constants. These are also ho- 
mogeneous universes with a 4-dimcnsional group acting 
transitively on the spacetime. Both the mean Hubble 
expansion rate and the shear are constant. Again, we 
see that the solution inflates anisotropically and is sup- 
ported by the existence of (3 / 0. It exists when a = 
and A = but not in the limit (3 — » 0. 



IV. AVOIDANCE OF THE NO-HAIR THEOREM 

The no-hair theorem for Einstein gravity states that 
for Bianchi types I — VIII the presence of a positive 



cosmological constant drives the late-time evolution to- 
wards the de Sitter spacetime. An exact statement of 
the theorem can be found in the original paper by Wald 
|| . It requires the matter sources (other than A) to 
obey the strong-energy condition. It has been shown 
that if this condition is relaxed then the cosmic no-hair 
theorem cannot be proved and counter-examples exist 
III LLJ- |33J, |3_^|. In [35(, the cosmic no-hair conjecture 
was discussed for Bianchi cosmologies with an axion field 
with a Lorentz Chern-Simons term. Interestingly, exact 
Bianchi type II solutions, similar to the ones found here, 
were found which avoided the cosmic no-hair theorem. 
However, unlike for our solutions, these violations were 
driven by an axion field whose energy-momentum tensor 
violated the strong and dominant energy condition. The 
no-hair theorem for spatially homogeneous solutions of 
Einstein gravity also requires the spatial 3-curvature to 
be non-positive. This condition ensures that universes 
do not recollapse before the A term dominates the dy- 
namics but it also excludes examples like that of the 
Kantowski-Sachs S 2 x S 1 universe which has an exact 
solution with A > which inflates in some directions 
but is static in others. These solutions, found by Weber 
[3(| . were used by Linde and Zelnikov [37j to model a 
higher-dimensional universe in which different numbers 
of dimensions inflate in different patches of the universe. 
However, it was subsequently shown that this behaviour, 
like the Weber solution, is unstable 38, 39] . We note that 
our new solutions to gravity theories with (3 ^ possess 
anisotropic inflationary behaviour without requiring that 
the spatial curvature is positive and are distinct from the 
Kantowski-Sachs phenomenon. 

The Bianchi type solutions given above inflate in the 
presence of a positive cosmological constant A. However, 
they arc neither de Sitter, nor asymptotically de Sitter; 
nor do they have initial singularities. Let us examine 
how these models evade the conclusions of the cosmic no- 
hair theorem. Specifically, consider the type II solution, 
ea. (|10fl We define the time-like vector n = d/dt orthog- 
onal to the Bianchi type II group orbits, and introduce 
an orthonormal frame. We define the expansion tensor 
@nv = n ii-,v and decompose it into the expansion scalar, 
9 = 9^ and the shear, a^ v = 9^ v - (l/3)(g M1/ + n M n„), 
in the standard way. The Hubble scalar is given by 
H = 9/3. For the type 7T metric, we find (in the or- 
thonormal frame) 

9 = 2b, = -diag(0,26,-6,-6). 

6 

As a measure of the anisotropy, we introduce dimension- 
less variables by normalizing with the expansion scalar: 



3a 



1 



1 1 



= diag(0, 



Interestingly, the expansion-normalised shear compo- 
nents are constants (and independent of the parameters 
a, (3, and A) and this shows that these solutions vio- 
late the cosmological no-hair theorem (which requires 
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G\ivl® — ► as t — > oo). To understand how this solu- 
tion avoids the no-hair theorem of, say, ref. [(J, rewrite 
eq.(0) as follows: 

In this form the higher-order curvature terms can be 
interpreted as matter terms contributing a fictitious 
energy-momentum tensor T^„. For the Bianchi II so- 
lution we find 

fpu = idiag (5b 2 - a 2 , -3b 2 + 3a 2 , -7b 2 - a 2 , -7b 2 - a 2 ) 
= diag(p,Pi,P2,P3)- (14) 

where p and pi are the energy density and the principal 
pressures, respectively. The no hair theorems require the 
dominant energy condition (DEC) and the strong energy 
condition (SEC) to hold. However, since p + pi + p2 + 
P3 = —3b 2 < the SEC is always violated when b ^ 0. 
The DEC is violated when p < and the weak energy 
condition (WEC) is also violated because p + p 2 = p + 
P3 = —{a 2 + b 2 )/2 < 0. These violations also ensure that 
the singularity theorems will not hold for these universes 
and they have no initial or final singularities. 

Are these solutions stable? Due to the complexity of 
the equations of motion it is difficult to extract informa- 
tion about the stability of these non-pcrturbative solu- 
tions in general. In the class of spatially homogeneous 
cosmologies the dynamical systems approach has been 
extremely powerful for determining asymptotic states of 
Bianchi models. A similar approach can be adopted to 
the class of models considered here; however, the com- 
plexity of the phase space increases dramatically due to 
the higher-derivative terms. Nonetheless, some stabil- 
ity results can be easily extracted. Consider, for exam- 
ple, a perfect fluid with a barotropic equation of state, 
p = wp, where w is constant. Due to the exponen- 
tial expansion, the value of the deceleration parameter 
is q = -(1 + H/H 2 ) = -1 for the type 77 and VI h 
solutions given. Hence, these vacuum solutions will be 
stable against the introduction of a perfect fluid with 
w > —1. This includes the important cases of dust 
(w = 0), radiation (w = 1/3) and inflationary stresses 
(-1 < w < -1/3). 

For perturbations of the shear and the curvature, the 
situation is far more complicated. Even within the class 
of Bianchi models in general relativity a full stability 
analysis is lacking. However, in some cases, some of the 
modes can be extracted. Consider again the Bianchi type 
II solution, eq. l(TTj|) . Using the trace of the evolution 



equations, eq.JSJl, we consider a perturbation of the Ricci 
scalar: 

R Ri 4A + r ie Xlt + r 2 e X2t . 

Using DR — — (R+8R), which is valid for spatially homo- 
geneous universes, eq.© again implies, to lowest order: 



ax, 2 --— ^ ± yi- 9F2(3a + /3) j, 

for (3a + /3) ^ 0. This shows that the perturbation of the 
Ricci scalar gives the same eigenmodes for the anisotropic 
solutions of type II and VIh as it did for perturbations 
of de Sitter spacetime in eq.(0. In order to determine 
the stability of other modes, like shear and anisotropic 
curvature modes, further analysis is required. 



V. DISCUSSION 

The solutions that we have found raise new questions 
about the thermodynamic interpretation of spacctimcs. 
We are accustomed to attaching an entropy to the geo- 
metric structure created by the presence of a cosmolog- 
ical constant, for example the event horizon of de Sitter 
spacetime. Do these anisotropically inflating solutions 
have a thermodynamic interpretation? If they are sta- 
ble they may be related to dissipative structures that 
appear in non-equilibrium thermodynamics and which 
have appeared been identified in situations where de Sit- 
ter metrics appear in the presence of stresses which vi- 
olate the strong energy condition 0, H2, H3, . They 
also provide a new perspective on the physical interpre- 
tation of higher-order gravity terms in the gravitational 
Lagrangian. 

In summary: we have found exact cosmological solu- 
tions of a gravitational theory that generalises Einstein's 
by the addition of quadratic curvature terms to the ac- 
tion. These solutions display the new phenomenon of 
anisotropic inflation when A > 0. They do not approach 
the de Sitter spacetime asymptotically and provide ex- 
amples of new outcomes for inflation that is driven by a 
p = —p stress and begins from 'general' initial conditions. 
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